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SUMMARY 


The  method  devised  by  W.  Vanderkulk  for  the  relative  location  of  three  hydrophones  in  “Remarks  on 
a  Hydrophone  Location  Method,”  JOURNAL  OF  UNDERWATER  ACOUSTICS,  Vol.  II,  No.  2,  April  1961 , 
yields  redundant  data  for  the  interhydrophone  distances  when  applied  to  an  entire  array.  The  purpose  of 
this  report  is  to  develop  the  mathematics  for  a  least-squares  method  to  use  the  redundant  measurements  to 
yield  the  coordinates  of  all  the  hydrophones  within  the  array  at  the  Pacific  Missile  Range  Facility,  Hawaiian 
Area  Underwater  Range. 
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INTRODUCTION 

The  method  developed  by  W.  Vanderkulk '  is  an  ingenious  technique  for  establishing  the  relative  loca¬ 
tion  of  three  hydrophones,  however,  problems  develop  when  an  attempt  is  made  to  extend  the  method  to 
locate  an  entire  airay  of  hydrophones. 

The  scheme  first  employed  was  to  calculate  the  interhydrophonr  distances  from  each  set  of  triads  and 
then  to  build  up  the  array  from  some  particular  starting  triad.  The  obvious  problem  is  that  since  the  results 
of  the  Vanderkulk  calculations  have  some  error,  the  common  distance  between  two  hydrophones  will  not  be 
the  same  when  measured  in  two  different  triads.  Thus  there  can  be  no  set  of  hydrophone  positions  that  will 
exactly  fit  all  the  measured  distances.  Some  arbitrary  judgment  must  be  made,  and  the  derived  hydrophone 
positions  must  be  constantly  revised  to  compensate  for  the  discrepancy  between  the  calculated  and  the  mea¬ 
sured  distances. 

As  would  be  expected,  the  errors  appear  to  be  more  pronounced  with  increasing  distance  from  the 
starting  triad. 

The  standard  technique  in  dealing  with  redundant  measurements  is  that  of  least  squares,  it  is  the  pur¬ 
pose  of  this  report  to  develop  a  least-squares  method  ot  using  the  Vanderkulk  information  to  simultaneously 
locate  the  entire  array  of  hydrophones. 


DEVELOPMENT  OF  THE  NORMAL  EQUATIONS 

The  first  step  in  a  least-squares  procedure  is  to  devise  the  function  to  be  minimized.  Since  the  most 
meaningful  output  from  the  Vanderkulk  calculations  is  the  interhydrophone  distances,  these  distances  will 
be  the  basis  for  the  minimization.  Thus,  let 

s=Iw,k,-D«)2  «> 

where  the  D|n[  represents  the  measured  distances,  the  represents  the  “true”  distance  between  pairs  of 
hydrophones,  and  the  W(  arc  the  1  elative  weights  of  the  i th  measurements. 


‘Department  of  Defense.  “Remarks  on  a  Hydrophone  Location  Method,”  by  \V.  Vanderkulk,  J.  UNDERWATER 
ACOUSTICS,  Vol.  II,  No.  2  (Apr  1961),  pp.  241-250. 


PRECEDING  PAGE  BUNK 
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The  next  step  is  to  derive  a  linear  relation  between  the  hydrophone  coordinates  and  the  distance.  The 
following  derivation  will  consider  only  the  X-Y  grid  coordinates  and  not  the  full  three-dimensional  location 
problem. 

Let  the  coordinates  of  the  Uth  hydrophone  (I  represents  the  row,  J  the  column)  be  X^,  Y-;  then 


Now,  let  Xf>  be  an  estimate  of  X~  and  Y|n*  an  estimate  of  Y-,  or 


x«  -  xr  *  v 


Y..  =  Y.(.n)  +  A..v 

ij  ij  ijY 


Introducing  the  notation 


■>S!&  - 

yJW] 

'-xg '? 

t(Yf-Yg>)2 

D(A)  = 

\sx)2  + 

(AijY-AkCY)2 

D^N)(X) 

-  (xr 

xg>) 

D(N)(Y) 

a(N)  _ 

D^X) 

ij.kfi 

n(N) 

*j,k2 

<©■ 

D(N)(y) 

we  have 

DiJ,k«  "[te)2  *  2D<N,(X)  (V  ■  4^x)  +  20<N,W(auy  -  AkSy)  *  DV)]2 

Now,  assume  that  second-order  differences  are  negligible  and  expand  the  terms  under  the  radical: 

Dij,k2  ~DijNkfi  +  aijjc£  (AijX  '  Ak£x)  +  bij,kfi  (AijY  "  Ak£y)  (H) 
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Substitution  of  (II)  into  (I)  yields: 


I  WiJ,k«{Dm(ij,k«)  -  [l®  *  ®  (v  - \tx).  *  ®  (4liV  -  \ty) j  (“0 


Thus,  taking  the  partial  derivatives  of  S  with  respect  to  the  A.jX  and  A~Y  and  equating  to  zero  will  yield 
the  standard  set  of  normal  equations. 

The  solution  to  the  normal  equations  will  be  the  best  estimate  of  the  hydrophone  X-Y  coordinates  in 
the  lm-’arized  least-squares  sense. 


DERIVATION  OF  THE  INITIAL  ESTIMATE 

The  system  of  equations  generated  by  the  minimization  of  equation  (III)  involves  an  iterative 
technique. 

The  derivation  that  follows  is  the  development  of  a  good  first  approximation  and  the  resulting  set  of 
equations  for  the  particular  configuration  at  the  Pacific  Missile  Range  Facility,  Hawaiian  Area  Underwater 
Range.  This  range  can  nominally  be  represented  as  a  collection  of  52  equilateral  triangles  with  sides  of 
8,000  feet  (see  figure  1).  The  center  of  the  range  is  hydrophone  34,  and  the  positive  Y-axis  is  determined 
by.  the  line  from  hydrophone  34  to  hydrophone  3-7. 


The  nominal  location  of  the  X-Y  coordinates  for  each  hydrophone  IJ  is 


(2r-l),s 

= 

(2-r)V3d,(s-4)d 

(2r,s)  = 

;  2  3Ms-  t)* 

(r  =  1 ,2,3)  (s  =  1  ,...,7) 
(r  =  1 ,2)  (s  =1 . 8) 


where 


d  =  8,000  feet 


and 


r  =  row 
s  =  column 

Now,  the  true  location  of  each  hydrophone  can  be  represented  by  this  nominal  location  and  some  additional 
displacement: 


(2r-l),sj  = 
(2r,s)  = 


(2  "  r)  +  e(2r-l)sx’(S  ‘  4)d  +  e(2r-l)sy 

^-rV3'd  +  e2rsx.(s--|)d  +  e2rsyJ 
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In  terms  of  these  coordinates,  the  distance  between  any  two  possible  adjacent  hydrophones  is 


D(2r-l)s,(2r-l)(s+l)  \/(e(2r-l)sx  '  e(2r-l)(s+l)x)2  +  "  e(2r-l)(s+l)y  +  e(2r-l)y^ 


^(2r-l)s,2rs  (  2  d  ’  e(2r-l)sx)  *(2  ~  e2rsy  +  e(2r-l)sy j 


^(2r-l)s,2r(s+l)  \j[  2  d  '  e2r(s+l)x  +  e(2r-l)sx)  +  (‘  2_e(2r-l)sy) 


D 


2rs,2r(s+l)  =  ^(e2rsx  '  e2r(s+l)x)2  +  f1  ‘  e 


2r(s+l)y  +  62rsy 


V3 

D2rs,(2r+l)(s-l)  =  \j[  2  d  '  e(2r+l)(s-l)x  +  e2rsx, 


e(2r+l)(s-l)y  +  e2rsy/ 


D 


2rs,(2r+l)s 


N 3 


9  d*e(2r+l)sx  +  e2rsy 


Hi 


e(2r+l)sy  +  62rsyj 


Next,  using  a  linearized  expansion, 


D/'i.  —  d-  (st').  1  " e, 


(2r-l)s,(2r-l)(s+l)-u'  \c(2r-l)sy  '  c(2r-l)(s+l)y) 
l(2r-l)s,2rs“d  +  2  f(2r-l)sx  ~e2rsx)  +  2  (f(2r-l)sy  ”  e2rsj^ 


D(2r-l)s,2r(s+ 1 )  d  +  ^2  f(2r-l)sx  ‘  e2i(s+l)y) 


D„„ d-  ^2rsy  “  e2r(s+l)y) 

vs , 


2rs,2r(s+l) 
)2rs,(2r+l)(s-l) 


iA-d+  ^  Ujjj ^  ^  ^ 


D 


2rs,(2r+l)s' 


V3/e 

2  \e2rsx’c(2r+l)(s-l)x/  ‘  2  c(2r+l)(s-l)y 
d  +  "^(€2rsx  ■e(2r+l)sx)  *^(f2rsy  "e(2r+l)sy) 


Thus, 


D|P)g  =  d  =  8,000  for  all 


V3  n 

a  =  ~y  or  0 


b  =  ±2 
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And  the  resulting  set  of  equations  (for  w;  =  1  for  all  weights)  is 


2  Ux 


V3, 


21x 


"^e22\~  2eHy  '  2£21y  +  e22y 


2^lsx -V^2sx  *N^e2(s+])x  'e2sy +  e2(s+l)y 

=  [Kis  ‘  d)  +  (D 1  s^2s  ‘  d)  +  (D11s?2(s+1))  +  (D1^2(s+1)  '  d)j 
(s  =  2,...,6) 

x/3^  r?  V3  ,1  ,  ,1 

3"2e17x"^e27x“T  e28x+2£17y  e27y +  2£2i  ■ 


(1) 


(2  through  6) 


"  [(^  -  d)  +  (D<17?27  -  d)  +  -  d)]  W 

~^T  e(2r-l)x  +^e2rlx  -^Te(2r+l)lx  e(2r-l'ly  H  e(2r+l)ly 

=  [ (D( 2r)  1 ,( 2r+ 1 ) l)  "  (°(*-l)l,(2r)l  “  d)J 

(r  =  1 ,2)  (8,21) 

"'ft e(2r-l)(s-l)x  e(2r-l)sx  +^'^3e2rsx  " ^  e(2r+l)sx  +  e(2r-l)(s-l)y  -e(2r-l)sy 
-e(2r+l)(s-l)y  +  e(2r+l)sy 


=  [(D2rs,(2r+l)(s-l)“d)  +  (D2rs,(2r+l)(s-l)  " d)  +  (D2rs,(2i+l)s  ~  d) 

+  (D2rs,(2r+l)s~d)  "  (D(2r-l)s,2rs  “  d)  “  (D(2r-l)s,2rs  “  d) 

~  (D(2r-l)(s-l)2rs  ~  d)  -  (D<2i-l)(s-l),2rs  “  d)J 
(r  =  1,2;  s  =  2,...,7)  (9  through  14  and  22  through  27) 
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2~  e(2r-l)7x  +  ^e2r8x  "  2  6(2r+l)7x  +  2  e(2r-l)7y  “  2  6(2r+l)7y 

=  (D2r8,(2r+1)7  _d)  “  (D(2r-l)7,2r8  “  d)j 

(r  =  1 ,2)  (15,  28) 

e21x  e22  x  +  3  V3  e31x  e41x  ->^e42x  +fe21y  ~e22y  "  2€41y  +e42y 

"  [(D3Ml-d)  +  fe-d)  +  -d)  -  fel-d)  -  fel-d) 


-  (<3,-4 

(16) 

-y/3  e2sx  -^e2(s+l)x  +4v/3"e3sx  -n^64sx  '^3e4(s+l)x+e2sy  _e2(s+l)y 

-e4sy  +  e4(s+l)y 

■  [(<4,-4  *  (<4,-4  *  (<4„4,)-4  +  (<4(,*l>-4 

-  (<3,-4 

-  (<3s"4  -  (D2(^!),3s  "  4  "  (D<3<s+„.3,-4] 

CO 

of 

H 

(17  through  20) 

\/T  nr  -  \/T  1  1 

e41x  -V3  e42x  +  3”e51x  +  2e41y  “e42y  +  2eSly 

=  [(<,s,-4*«s,-4  +  (<.»-d)] 

(29) 

“■^S4sx  ~^3  e4(s+l)x  +  2N^e5sx  +e4sy  _e4(s+l)y 

■  [-(<4|,,Ss-4  -  (D4(s+1),5s  -  4  *  (<5,-4  *  (<5,-d 

')l 

(s  =  2,. ..,6) 

(30  througli  34) 

-  \/3  e47x  -  3  e48x  +  3  e57x  +  e47y  “  I  e48y  ‘  K?y 

"  (D48>,57  ~  d)  “  (D47),S7“dj  “  (Dv?\si  “  d)j  (35) 
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Z5-.V3.V3  .7  ,  1 

■V3  ellx  "  4  e21x+  2  e22x  +  4elly  '  fc12y  "e21y  '  2  e22y 

■  -Hi!.2  -  <*) +  #'&.  -  o)  -  - 1)  -  RH)]|  (36) 

Vf  \/3  2  1 

2  e2«  2  e2(s+l)x‘el(s-l)y +  3elsy 'el(s+l)y  ‘  262sy  ‘  262(s+l)y 

=  IRmu.  -  <■) '  R’lcm) ' d)  '  jIKW.,  - d) +  - d)  -  Ri,  - ») 

(s  -  2,. ..,6)  (-37  through  41) 

V?  \/3  V3  7  2  - 

4  e17x  4  e27x  2  628x  '  e16y  +  4  e17y  '  2  f27y  '  4  e28y 

■  (-Rn"*)  4|te8-d)-R2,-d)-te,-d)]|  (42) 

V3  V3  2  5  1 

4  e(2r-l)lx  4  e(2r+l)lx  '  4  e(2r-l)ly  +  2  e(2r)ly '2  e2r2y  "  4  e(2r+l)ly 

=  j-(D(2r)l,(2r)2  -d)  '  (D(Ir)l  ,2r2  *  d)  +  j[(D(2r-l)l  ,(2rl) ' d)  +  (D(2r)  1  ,(2rH)  " d)]  j 

(43,  57) 


(r=l,2) 


V3 


.  VI 


V3 


1 


Y  e(2r-l)(s-l)x  ‘  ~  e(2r-l)sx  '  “  e(2r+ 1  )(s-l )x  +  ~  £(2r+l)sx  '  2  e(2r-l)(s-l)y 
'  2  6(2r-l)sy  '  2  e2r(s-l)y  +  6  e(2r)(2+l)y  '  2  e(2r+l)(s-l)y  *  2  e(2r+l)sy 

=  |  (D2r(s-l),2rs  '  d)  +  (D2r(d-1  ),2rs  *  d)  *  (D2rs,2i(s+ 1) ' d)  *  (D2rs,2r(s+ 1 )  *  d) 

‘  2l(D(2r-l)s,2rs  '  d)  +  (D(2M)s,2rs  "  d)  +  (D2rs,(2r+l)s  '  d)  +  (D2rs,(2r+l)s  '  d) 
*  (D(2r-l)(s-l),2rs  '  d)  '  (D(2r-l)(s-l),2rs  *  d)  *  (D2rs,(2r+l)(s-l) ' d) 
'(D2rs,(2r+l)(s-l)*d) 

(r  ~  1  >2;  s  =  2,.. .,7) 


(44  through  49  and  58  through  63) 


\/T  %/3  1  5  l 

4  e(2r-l)7x‘  4  6(2r+l)7x  '  4e(2r-l)7y e2r7y  +  2€2(4+l)y  '  4e(2r+l)7y 

"  {(D2r7,2r8  ‘  d)  +  (D2r7,2r8  '  d)  +  2[(D(2M)7,2r8  '  d)  +  (D2r8,(2r+1)7  '  d)  } 


(r=l,2) 


(50, 64) 


vll,  J\/3,  vfc  ,VIC  i  1  ,  7  .  1  1 

4  21x  4  e22X  4  e41x  2  e42x  '  4e21y  "  2  e22y  2  e31y 632y  '  4  e41y  '  2C42y 

=  I'fe  -  «)  -  K’,32  -  i)  -  K’,„  -  d)  ♦  (Dg>,„  -  d)  ♦  (D<  V,„  ■  d)  *  (d®2  -  d) 

”  (D3M1  ’d)  “  (D2K31  '^)j  <51 


V3  \/3  \/3  V3  1  1 

2  e2sx  "  2  £2(s+1)x  “  2  £4sx  +  “  e4(s+l)x  "  2  e2sy  *  2  e2(s+l)y  ‘^e3(s-l)y  +  de3sy 


1  1 


^  * 
'ze3(s+l)y  ‘  264(s+l)y 


=  1(D3(s-1),3s  "  d)  +  (D3(s-1),3s  *  d)  *  (^(s+l) ’ d)  '  (D3s?3(s-1) ’ d)  "  4 [(D2(s+ 1 ) ,3s  *  d) 
+  (D2(s+1),3s  '  d)  +  (D3s^4(s+ 1) ' d)  +  (D3$!4(s+1)  *  d)  *  (D2s!3s  *  d)  ‘  (D2s^3s  '  d) 
-K!4,-d)- (D<3)4s-d)H 

(s  =  2,3,5, 6)  (52  through  55) 

n/3  VI  \/3  V?  2  1  7  11 

2  e27x  ‘  4  e28x  '  2  647x  +  4  e48x  '  2  e27y  ’  4  e28y  "  2  e36y  +  2  e37y  '2  e47y  "4  e48y 

=  ' d)  +  (D36.37  '  d)  '  2  '  d)  +  ' d)  '  ”  d)  '  (D27),37  *  d) 

'  (D3V,47  '  d)  ‘  {^37^,47  *  d)| }  (56) 

V3  c  V3  V3  2  1  .  7  _ 

4  e41x  "  2  42x  4  e51x"4e41y'2e42y  4e51y'2€52y 


=  |  '(D51),52  '  d)  *  2  [  (D42^,5 1  ‘  d)  +  (D42,51  *  d)  ’  (D41},51  *  d)]  | 


n 


y/3  .  s/3  l  l  , 

2  e45x  ■  2  e4(s+l)x  "  2  e4sy  '  2  e4(s+l)y  "  e5(s-l)y  +  i  e5sy  ~e5(s+l)y 

=  I  (DsU-D.Ss  -  d)  -  1 }  -  d)  -  i  i  (DV(«  1 )  ,5s  -  d)  +  (D4<U  1  ).5s  -  d) 


(s  =  2,... ,6) 


1 


_  _1  .7 

2  e47x  ”  4  e48x  '  4  e57x  '  2  647y  *  4  e48y  ‘  e56y  4  e57y 


\/5  ,  V3  V3 

l-7x  ”  4  e48x  *  4 

-  [(“$«  -  »)  -  j(D«!57  '  »)  -  K!s7  -  d)  -  (»$57  '  »)] 


(66  through  70) 


(71) 


Thus, 


X(D  =  x(0)  +  e. 


1JX 


yd)  =  Y(0>  +  e. 

ijy 


and 


and  the  iteration  is  continued  until  |  -  XKN+1^  |  <  e, .  |  -  y[jN+1^  [  <  e2  for  all  i,j. 

Px 

To  further  simplify  the  program  it  was  also  assumed  that  aj^r^  =  constant  =  ±—  or  0,  and 
=  constant  =  or  0.  The  quality  of  the  results  obtained  imply  that  this  assumption  is  valid,  and 

1J>KX  2 

it  does  not  seem  necessary  or  worthwhile  to  include  the  non-constant  portion  of  the  coefficients  in 
the  calculation. 

The  system  may  be  written  more  concisely  in  matrix  notation 

A  /a  A  \ 

AAX  =  (b-b(N)j 

where 

A  =  the  matrix  of  coefficients 
A 

AX  =  the  vector  of  displacement  from  the  previous  best  estimate  of  position 
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A  A 

b  =  a  vector  composed  of  linear  combinations  of  the  interhydrophone  distances  (d) 


b(N)|=  the  vector  of  corresponding  Nf/i  estimates 

A  further  consequence  of  the  linear  least-squares  formulation  of  the  problem  is  an  estimate  of  the 


propagated  error  in  the  coordinates  X  =  X(N  -  1)  +  A-1  |b  -  b(N  -  1)J .  Thus  it  can  be  shown  (see  appendix) 
that  oA  =  A-1  ^A_1j  .  Further,  b  =  Cd,  so  that,  finally,  aA  =  A-1C  oA  ^A_1cj  ,  and  for  an  estimate  of 

aA  =  DIAG  j ,  J56 j ,  taking  the  diagonal  elements  of  will  give  an  estimate  of  propagated 

errors  in  each  of  the  hydrophone  coordinates  resulting  from  the  uncertainty  in  the  interhydrophone  dis¬ 
tance  measurements. 


SUMMARY  OF  RESULTS 

Although  only  one  Vanderkulk  Survey  has  been  run  at  the  Pacific  Missile  Range  Facility,  Hawaiian  Area 
Underwater  Range,  four  sets  of  interhydrophone  distances  were  used  as  input,  resulting  in  four  sets  of  hydro¬ 
phone  positions.  The  first  set  is  the  original  unedited  results  obtained  by  ITT  (International  Telephone  and 
Telegraph  Corporation).  The  second  set  is  the  ITT  data  with  certain  distances  edited.  The  ITT  results  with 
this  second  set  of  interhydrophone  distances  are  currently  in  use  at  this  range.  The  third  set  was  obtained  at 
the  Pacific  Missile  Range  by  the  Mathematical  Analysis  Branch  (Code  3433),  the  interhydrophone  distances 
were  recomputed  with  a  greater  amount  of  data  per  triad,  and  the  resulting  hydrophone  positions  were  recal¬ 
culated.  The  fourth  set  of  input  resulted  from  an  editing  of  the  Pacific  Missile  Range  data. 

With  these  four  sets  as  input,  four  least-squares  solutions  were  obtained.  Table  1  is  a  tabular  listing  of 
the  difference  between  the  results  of  the  least-squares  run  4  (taken  as  the  standard)  and  the  other  runs  for 
corresponding  coordinates.  The  difference  between  the  present  set  of  coordinates  and  the  proposed  set  is 
presented  graphically  in  figure  2.  As  would  be  expected,  the  results  from  the  least-squares  formulation  are 
much  more  consistent. 

Table  2  lists  the  estimate  of  propagated  error  in  each  coordinate  caused  by  the  uncertainties  in  the 
interhydrophone  distance  measurements.  The  average  value  is  about  5  feet. 

As  a  further  consistency  check,  the  interhydrophone  distances  were  computed  with  the  use  of 
the  calculated  hydrophone  positions.  Here,  the  results  are  quite  revealing;  the  difference  between  the 
interhydrophone  distances  computed  and  the  weighted  average  of  the  measured  distances 


d=  [ff2/(a1+^]d1  +  [®1 

/(ai+02 

jjd2  is  tabulated  below: 

Difference  (Feet) 

ITT<  1) 

LS(1) 

ITT(2> 

LS(2) 

PMR 

LS<3) 

<2 

60 

68 

66 

83 

63 

84 

2  to  5 

9 

19 

14 

5 

17 

4 

5  to  10 

6 

1 

8 

8 

10  to  15 

7 

15  to  20 

3 

Over  20 

3 

NOTE:  LS  =  least  squares. 
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Table  1.  Difference  in  Coordinate  Position  From  Those  of  ST 


Hydro¬ 

phone 

ITT(  1)-ST 

ITT(2)-ST 

PMR-ST 

SLSS(1)-ST 

SLSS(2)-ST 

SLSS(3)-ST 

11 

-30.65 

24.84 

5.86 

-4.39 

17.20 

-1..42 

32.71 

-0.32 

3.00 

7.26 

-0.18 

-0.30 

12 

-15.85 

20.40 

1.62 

-2.27 

11.38 

-1.26 

28.00 

-4.58 

3.80 

2.04 

-5.32 

-0.44 

13 

-8.44 

9.02 

-2.82 

-7.29 

0.39 

-2.00 

26.26 

-6.59 

3.99 

1.66 

-4.81 

-0.89 

14 

8.64 

7.97 

7.49 

-4.32 

5.47 

-0.48 

27.84 

-4.45 

5.16 

1.31 

-3.25 

-0.85 

15 

18.80 

11.18 

5.26 

-3.70 

7.46 

-0.96 

29.19 

-3.02 

4.75 

2.98 

-4.22 

-0.80 

16 

25.31 

6.55 

6.40 

-2.23 

3.90 

-1.46 

29.44 

-2.56 

4.84 

4.12 

-4.01 

-0.69 

17 

12.75 

-5.34 

-2.88 

-5.95 

-5.54 

-1.29 

34.66 

2.66 

5.83 

8.83 

0.47 

-0.45 

21 

-34.80 

29.11 

11.63 

-5.60 

19.06 

-1.62 

17.00 

0.03 

12.54 

-1.10 

-2.48 

-0.36 

22 

-21.32 

24.38 

3.82 

-2.60 

15.32 

-1.35 

15.41 

-1.30 

5.16 

3.22 

1.98 

-0.29 

23 

-13.24 

13.53 

-2.88 

-0.51 

9.82 

-1.26 

14.55 

-2.41 

3.38 

4.39 

1.45 

-0.15 

24 

-0.92 

7.35 

3.74 

-2.69 

3.90 

-0.10 

11.68 

-4.77 

0.48 

2.36 

-3.20 

-0.11 

25 

8.80 

10.60 

6.27 

-3.79 

9.21 

-0.69 

9.94 

-6.22 

5.17 

0.15 

-5.50 

-0.41 

26 

16.47 

6.67 

8.54 

-2.43 

6.29 

-1.22 

10.44 

-5.75 

5.04 

1.56 

-2.80 

-0.39 

27 

9.43 

0.99 

5.31 

-1.39 

1.81 

-1.56 

10.77 

-5.41 

5.00 

1.24 

-2.77 

-0.67 

28 

2.62 

-6.62 

1.53 

-1.31 

-3.74 

-0.82 

11.01 

-5.16 

3.26 

1.57 

-2.23 

-0.79 

31 

-34.04 

28.87 

7.95 

-4.22 

16.86 

-1.65 

-5.05 

4.74 

7.71 

1.23 

6.92 

-0.61 

32 

-20.57 

22.17 

1.83 

-1.26 

14.05 

-1.23 

-1.33 

8.58 

9.68 

2.47 

9.13 

-0.50 

33 

-5.28 

13.27 

1.17 

-0.10 

10.03 

-0.85 

1.60 

0.60 

-3.25 

2.30 

1.01 

-0.46 

34 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

35 

12.81 

12.70 

4.81 

-4.11 

6.90 

-1.26 

1.39 

1.40 

0.48 

0.32 

0.78 

-0.02 

36 

2.66 

3.87 

0.96 

-1.92 

3.96 

-1.13 

1.30 

1.32 

-0.42 

-0.06 

0.78 

-0.03 

37 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

2.14 

2.15 

-1.40 

0.57 

1.50 

0.16 

41 

-49.39 

16.03 

4.94 

-19.25 

2.96 

-1.58 

-2.24 

12.38 

0.24 

3.06 

11.66 

-0.57 

42 

-21.94 

24.80 

3.93 

-5.08 

13.21 

-1.15 

-4.24 

11.57 

4.01 

-1.21 

7.80 

-0.62 

43 

-8.86 

20.36 

1.72 

-0.51 

14.07 

-1.01 

-6.78 

8.98 

-0.59 

-2.72 

5.49 

-0.80 

44 

-2.95 

6.62 

-0.46 

-1.92 

5.67 

-0.83 

-6.49 

9.29 

-1.13 

-3.04 

5.75 

-1.01 

45 

2.89 

12.50 

0.43 

-2.22 

13.02 

-2.10 

-8.27 

7.51 

-1.75 

-2.05 

6.83 

-1.10 

46 

-6.01 

3.64 

-1.18 

-3.81 

4.74 

-1.27 

-9.51 

6.27 

-2.23 

-1.12 

5.68 

-0.70 

47 

-6.99 

2.69 

-0.09 

-0.89 

2.09 

-0.60 

-9.62 

6.16 

-2.07 

-1.29 

5.03 

-0.70 

48 

-13.75 

-4.03 

-3.45 

0.43 

-1.62 

0.24 

-9.58 

6.20 

-2.24 

-1.68 

5.11 

-0.73 

51 

-25.43 

21.02 

-3.36 

-11.65 

8.47 

-1.34 

-6.84 

8.42 

-0.63 

-8.74 

3.35 

-0.92 

52 

-7.18 

21.99 

2.91 

-2.36 

14.20 

-1.02 

-6.08 

9.57 

0.03 

-7.04 

4.80 

-0.87 

S3 

-0.57 

8.98 

-0.12 

-1.63 

9.18 

-0.86 

-3.85 

11.90 

-0.69 

-2.58 

8.90 

-0.87 

54 

-3.57 

6.02 

-0.07 

-6.51 

4.75 

-1.40 

-4.87 

10.88 

-1.57 

-0.06 

8.16 

-0.66 

55 

-0.61 

9,01 

0.14 

-3.44 

9.62 

-1.84 

-2.92 

12.82 

-1.43 

0.20 

11.16 

-1.17 

56 

-6.53 

3.13 

-0.96 

-2.97 

3.56 

-1.01 

-4.42 

11.32 

-1.83 

-1.13 

10.22 

-1.29 

57 

-9.55 

0.15 

-1.25 

0.34 

0.73 

-0.20 

-4.82 

10.93 

-0.67 

-2.50 

8.73 

-1.36 

Note:  SLSS(l)  =  Least-squares  results  from  same  input  as  used  for  ITT(l). 

SLSS(2)  =  Least-squares  results  from  same  input  as  used  for  ITT(2). 

SLSS(3)  =  Least-squares- results  from  same  input  as  used  for  PMR. 

ST  =  Least-squares  results  from  edited  PMR  input. 
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•  NOMINAL  LOCATION 
O  PRESENT  LOCATION 

*  PROPOSED  LOCATION 

Figure  2.  Comparative  Hydrophone  Coordinates. 
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Table  2.  Estimate  of  Propagated  Error 


Hydro¬ 

phone 

1  a 

Error  (Feet) 

Hydro¬ 

phone 

1  c 

Erro-  (Feet) 

11Y 

7.78 

33X 

3.64 

11X 

9.73 

35Y 

2.88 

12Y 

5.38 

35X 

2.88 

12X 

6.46 

36Y 

3.06 

13Y 

5.00 

36X 

3.05 

13X 

4.73 

37Y 

3.87 

KY 

4.67 

41Y 

4.68 

14X 

2.70 

4 IX 

11.63 

15Y 

4.30 

42Y 

4.09 

15X 

3.27 

42X 

7.51 

16Y 

4.54 

43Y 

3.69 

16X 

4.16 

43X 

4.95 

17Y 

4.79 

44Y 

3.31 

17X 

5.14 

44X 

2.80 

21Y 

5.91 

45Y 

3.27 

2 IX 

12.47 

45X 

3.28 

22Y 

3.73 

46Y 

3.31 

22X 

7.46 

46X 

3.22 

23Y 

3.26 

47Y 

3.40 

23X 

4.93 

47X 

2. 16 

24Y 

3.03 

48Y 

3.76 

24X 

2.51 

48X 

3.08 

25Y 

2.85 

S1Y 

6.96 

25X 

2.30 

SIX 

9.62 

26Y 

3.13 

52Y 

5.87 

26X 

3.35 

52X 

6.40 

27  Y 

3.66 

53Y 

5.58 

27  X 

4.37 

53X 

4.19 

28Y 

4.40 

54Y 

5.18 

28X 

4.66 

54X 

3.41 

31Y 

3.71 

55Y 

4.41 

31X 

3.95 

55X 

3.67 

32Y 

2.83 

56Y 

4.31 

32X 

6.06 

56X 

3.31 

33Y 

2.45 

57  Y 

5.13 

57  X 

2.75 

CONCLUSIONS 

A  study  of  table  1  and  figure  2  reveals  two  very  important  phenomena.  The  first  is  that  the  least- 
squares  solution  obtained  with  the  unedited  ITT  data  actually  agrees  with  the  assumed  standard  much  better 
than  the  least-squares  solution  obtained  from  the  edited  ITT  data.  The  second  is  that,  relative  to  the  assumed 
standard,  the  range  appears  to  be  slightly  rotated,  or  the  fourth  and  fifth  rows  appear  to  be  biased  in  the 
Y-direction;  the  fourth  row  by  about  8  feet  and  the  fifth  row  by  about  10  feet. 

As  was  previously  stated,  these  results  were  obtained  with  the  use  of  simplifying  assumptions,  and 
future  studies  are  planned  to  show  the  effect  of  different  weighting  schemes,  and  possibly  also  of  the  intro¬ 
duction  of  the  complete  coefficient  matrix  (including  the  functional  dependence). 

The  present  program  was  developed  to  test  the  feasibility  of  the  technique,  and  this  has  been  amply 
verified.  The  convergence  of  the  solution  is  such  that  after  five  iterations,  the  hydrophone  coordinates  are 
consistent  to  within  a  maximum  of  0.1  foot. 
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A  further  benefit  is  that  the  entire  run  (five  iterations)  is  completed  in  less  than  2  minutes  of  computer 
time  rather  than  the  days  or  weeks  required  by  the  present  method  that  includes  hand  calculations.  Also,  as 
previously  stated,  there  is  an  automatic  error  estimate. 

A  further  advantage  of  this  technique  is  that,  as  future  surveys  are  performed,  these  data  can  be  com¬ 
bined  with  previous  survey  information  to  yield  an  updated  best  estimate  that  would  include  all  available 
information. 

Obviously,  this  technique  can  be  applied  to  any  other  underwater  range.  The  only  difference  would 
be  in  the  formulation  of  the  initial  estimates. 

This  report  has  been  devoted  to  finding  only  the  X-Y  grid  coordinates  of  the  hydrophones.  Since  the 
Vanderkulk  survey  measures  the  depth  of  each  hydrophone  directly,  all  that  is  required  for  the  Z-coordinates 
is  a  weighted  average  of  these  measuremt  nts. 


17 


APPENDIX 

DERIVATION  OF  PROPAGATED  ERROR  ESTIMATE 


At  any  stage  of  the  iteration,  X(N)  -  X(N  -  1)  =  A'1  (b  -&(N  - 1))  ,  or,  in  the  component  form, 
Xj(N)-X,(N-I)=SA:>(bk-bk(N-l)). 


Thus, 


3x/N)  8Xj(N  -  1)  y  3bk(N-l) 

3bm  9bm  jk  *  ^  ^  abm 


Further, 


I 


a.i5s2L1> 

*  9bm 


I 


3bk(N  - 1)  9xe\ 
9xfi 


Now,  if  bk(x)  is  expanded  in  a  Taylor  series, 


bk(x)  =  bk(N  -  1)  + 


or 


9bk 

Aki="3xJ 


PRECEDING  PAGE  BLANK 
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Thus,  since 


0xj  _  0xj(N  “ !)  _j  3xj(N-l) 

0b^  -  9bm  +  Ajkm  3b~  =  Ajkm 


Therefore 


oa-A~1  o^A^y 


NOTE:  The  preceding  derivation  holds  only  when  A  is  assumed  to  be  a  matrix  of  constants;  when  the  full 
coefficient  is  used  a  slightly  different  formula  arises.  The  complete  formulation  appears  in  a  Pacific 
Missile  Range  unpublished  Working  Note  3433-21 ,  Propagation  of  Error  Associated  with  an  Iterative 
Least-Squares  Solution,  by  S.  Berman,  dated  December  1967. 
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